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ABSTRACT 

The emission of FU Orionis objects in outburst has been identified as arising in 
rapidly accreting protoplanetary disks, based on a number of observational properties. 
A fundamental test of the accretion disk scenario is that the differentially rotating disk 
spectrum should produce a variation of rotational velocity with the wavelength of ob- 
servation, as spectra taken at longer wavelengths probe outer, more slowly rotating disk 
regions. Previous observations of FU Ori have shown smaller rotation at near-infrared 
(~ 2.2/im) wavelengths than observed at optical (~ 0.6/im) wavelengths consistent with 
the assumption of Keplerian rotation. Here we report a spectrum from the Phoenix in- 
strument on Gemini South which shows that differential (slower) rotation continues to 
be observed out to ~ 5/xm. The observed spectrum is well matched by the prediction 
of our accretion disk model previously constructed to match the observed spectral en- 
ergy distribution and the differential rotation at wavelengths < 2.2fim. This kinematic 
result allows us to confirm our previous inference of a large outer radius (~1 AU) for 
the rapidly accreting region of the FU Ori disk, which presents difficulties for outburst 
models relying purely on thermal instability. While some optical spectra have been 
interpreted to pose problems for the disk interpretation of FU Ori, we show that the 
adjustment of the maximum effective temperature of the disk model, proposed in a 
previous paper, greatly reduces these difficulties. 

Subject headings: accretion, accretion disks — stars: formation — stars: pre-main 
sequence 



1. Introduction 

FU Orionis systems are a class of exceptionally luminous young stellar objects found in star- 
forming regions (Hartmann & Kenyon 1996). Originally identified from their very large increases in 
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optical brightness over timescales of years or less (Herbig 1977), a larger group of heavily extincted 
probable members of the class have been identified by their characteristic infrared spectra, which 
strongly differ from those of typical T Tauri stars (Reipurth &: Aspin 1997; Aspin & Reipurth 
2003; Reipurth et al. 2007). Additional support for the identification of these heavily embedded 
objects comes from recent high-resolution infrared spectroscopy, which indicates that many of these 
sources are rapidly rotating and exhibit double-peaked absorption line profiles as observed in FU 
Ori (Greene, Aspin, & Reipurth 2008). 

The accretion disk model for FU Ori objects (Hartmann & Kenyon 1996, and references therein) 
rests fundamentally on the need to explain the peculiar spectral energy distributions (SEDs) of these 
objects, which are much broader than that of a single temperature blackbody or star, and which 
exhibit a continuously varying spectral type as a function of wavelength. The disk model naturally 
accounts for these properties, as observations at longer wavelengths probe increasingly cooler disk 
regions with later spectral types; our detailed model for FU Ori matches the SED from optical 
wavelengths to the mid-infrared region (Zhu et al. 2007, 2008). In addition, the disk model predicts 
that differential rotation should be observed, with slower rotation seen at longer wavelengths, which 
arise from outer disk radii, and this has been confirmed by comparing optical 0.6^m) and near- 
infrared (~ 2.2/im) spectral line profiles (Hartmann & Kenyon 1987a, b; Kenyon, Hartmann, & 
Hewett 1988). 

A desirable feature in a model or theory is an ability to make predictions that can be tested 
observationally. The disk model predicts that the observed rotational spectral line broadening 
should be even smaller at A ~ 5//m (using the fundamental CO vibrational transitions) than at 
2.2/um (using the first overtone CO vibrational lines). In this paper we present a high-resolution 
spectrum of FU Ori in the 5//m region which matches the predictions of the disk model. We also 
show that some discrepancies seen in optical spectra of FU Ori in comparison with simple disk 
models are alleviated by a decrease in the maximum disk temperature which we proposed in Zhu 
et al. (2007). 



2. Observations 

A high-resolution spectrum of FU Ori at 4.9/um was obtained at UT 01:00:29 on 2007 February 4 
using the Phoenix spectrometer (Hinkle et al. 1998, 2000, 2003) on the 8-m Gemini South telescope. 
Observations were taken with a two-pixel slit (0.17") for a resolution of X/5X = 75,000 over the 
wavelength range 4.808 — 5.050^m. We observed FU Ori at two positions along the slit for eight 
2 minute exposures. We also observed the B2 III star HR1790 for telluric line correction and took 
10 flat-field and dark images. 

We reduced the data using IRAF.0 We averaged the flat-field and dark images and subtracted 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 



- 3 - 



the average dark image from the average flat-field image. This averaged, dark-subtracted flat- 
field image was then divided into the target spectra. Images at different positions of the slit were 
differenced to remove the sky and dark backgrounds. We then extracted the spectra using the IRAF 
apall routine and later combined and flattened the spectra using splot in IRAF. The spectrum of 
the hot star HR1790 was used to divide out the telluric lines from the FU Ori spectrum. Wavelength 
calibration was computed using telluric lines from the Arcturus atlas of Hinkle et al. (1995). 

We also obtained a high-resolution spectrum of FU Ori at UT 17:59:15 on 2008 November 
14 using the Magellan Inamori Kyocera Echelle (MIKE) spectrograph on the 6.5 m Magellan 
Clay telescope at Las Campanas observatory (Bernstein et al. 2003). MIKE is a double echelle 
spectrograph which delivers full wavelength coverage from about 3350 — 5000 A (blue side) and 
4900 - 9500 A (red side). The data were obtained in subarcsecond seeing with a 0.7"slit, binning 
2x2 and an exposure time of 120s. The resolutions were ~ 40,000 and ~ 30,000 for the blue and 
red sides respectively. The MIKE data were reduced using the MIKE Redux IDL pipeline. @ 



3. Results 

Figure 1 shows the reduced Phoenix spectrum of FU Ori. Strong telluric features mean that 
regions around wavenumber 2009.3 cm~^ are not usable, and small residuals from the telluric 
correction can be seen at wavenumbers 2010.4, 2010.9, 2011.9, and 2013.5 cm^^. The spike around 
wavenumber 2012.5 cm~^ is due to a bad pixel. Outside of these regions, an absorption spectrum 
is clearly present, with relatively broad features and substantial blending. Most of the lines are due 
to the P-branch of the fundamental rotational- vibrational transitions of CO. 

To interpret the results, we calculated a synthetic disk spectrum using the methods described 
by Zhu et al. (2007). The model parameters were those used by Zhu et al. (2007) to fit the 
SED of FU Ori and to match the observed rotational broadening observed at optical and near- 
infrared wavelengths: central star mass ~ 0.3 Mq, mass accretion rate of the inner high M disk 
~2.4x10~^Mq yr~^, disk inner radius ~5 Rq, outer radius of the inner high M disk ~ 1 AU, and 
the disk inclination angle ~55°. (See Figure 8 of Zhu et al. for the fit to the 2.2/im CO lines). We 
emphasize that we have not changed or adjusted any parameters from the Zhu et al. (2007) FU Ori 
model; these are predicted spectra. 

The lower dotted curve in Figure 1 shows the synthetic disk spectrum observed pole-on, so that 
individual spectral lines can be seen without the blending that occurs due to the large rotational 
broadening. Comparison of the nonrotating spectrum with the observations shows that FU Ori 
has significantly larger line widths, consistent with rapid rotation, and unlike profiles of M giants 
and supergiants. The upper dotted curve shows the synthesized spectrum using the inclination and 
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central mass used to obtain a match to the 2.2fj,m CO Hne widths. The agreement between synthetic 
and real spectra is quite good, except near 2010. 7cm~^ where we are missing the CO 5-4 P7 line in 
the modd: Th. could be due to ad;pting too .majl an oscillator strengthl The half-width at half- 
depth (HWHD) of the lines in this spectral region (5 /im) is ~ 22kms~^, considerably smaller than 
the line widths measured at 2.2/im (HWHD ~ 36 km s~^) (Hartmann &; Kenyon 1987a; Hartmann, 
Hinkle, & Calvet 2004; Zhu et al. 2007). This HWHD at 5 /im is close to the Keplerian velocity at 
0.5 AU around 0.3 Mq central star. All the strong lines have been identified in the pole-on model 
spectrum. Though some water lines are present in this spectrum, they are washed out or blended 
in the broadened spectrum (upper dotted curve). Only the strong CO fundamental lines can be 
identified with certainty. 

It is also worth noting that some of the unblended line profiles exhibit evidence for double- 
peaked shapes predicted by simple disk models (e.g., Hartmann & Kenyon 1987a,b). In addition, 
some of the blends show sharp features (eg. 2011.2 cm~^ line) which are the result of overlapping 
double-peaked lines. These features naturally arise in a disk model but would not be seen in 
rotating star models (unless large polar spots are invoked; see below). 

In Figure 2 we display a segment of the MIKE spectrum of FU Ori in the wavelength range 
7030 — 7100 A for comparison with the synthetic Keplerian disk spectrum. We again find good 
agreement between model and observation, demonstrating that there has been no change in the 
estimated optical rotational velocity of the object between the observations in Zhu et al. 2007, 
which we used to set the disk parameters, and this paper. The HWHD of the optical lines in this 
wavelength range is ~65ib5 km s~^, consistent with HWHD ~ 62ib5 km s~^ measured by Petrov 
& Herbig (2008; PH08). 

Thus, compared with HWHD of 2 micron CO first-overtone lines ^ 36±3 km s~^ and HWHD 
of 5 micron CO fundamental lines ~ 22zb2 km s~^, the differential rotation in FU Ori observed 
over nearly an order of magnitude in wavelength is consistent with Keplerian rotation. The slow 
rotation observed at 5//m implies spectral formation at radii out to 0.5 AU, in agreement with 
our SED modeling; this supports our conclusion in Zhu et al. (2007, 2008) that the extent of the 
high-accre tion rate disk is larger than can be explained by pure thermal instability models for 
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4. Discussion 

The consistency of the variation of rotational velocities as observed over to A ~ 0.7 — 5//m 
with Keplerian rotation seemingly provides strong evidence for the accretion disk interpretation of 
FU Ori. However, PH08 recently argued that while the infrared spectrum is that of an accretion 
disk, the optical spectrum is produced by a rapidly-rotating star with a dark polar spot. The 
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PH08 argument against a pure disk model for a central star rests on three main points: there is 

no evidence for a variation of absorption line width as a function of excitation potential over the 
wavelength range A ~ 0.52 — 0.86;um; there is no evidence for a variation of rotational velocity 
with wavelength over that wavelength range; and the observed line profiles are more "boxy" or 
fiat-bottomed than the double-peaked profiles of the disk model. 

It is important to recognize that the above-listed effects expected for a disk spectrum require 
not only differential rotation but a temperature gradient as well. A Keplerian disk exhibiting a 
constant effective temperature would not show any effect of rotational velocity with either excitation 
or wavelength; and the double-peaked behavior of line profiles only occurs because the outer, 
slowly rotating regions do not fill in the profile at line center, as these regions are too cool to 
emit significantly at the wavelength of observation. While the standard steady disk temperature 
distribution T^ff oc [1 — is a power law at large radii, it is relatively flat at distances 

within about twice the inner radius Ri. Therefore, observations at long wavelengths which probe the 
outer disk where the temperature falls rapidly with radius will exhibit stronger rotational velocity 
variations and more double-peaked line profiles than observations at short wavelengths probing the 
inner, more nearly isothermal disk. 

In Zhu et al. (2007) we were forced to use a maximum disk temperature T^aa; = 6420 K to 
match the SED of FU Ori, which is lower than the 7200 K maximum temperature adopted in the 
model used by PH08 (which was based on the earlier model by Kenyon, Hartmann, & Hewett 1988). 
Lowering the maximum temperature has the effect of making the flatter part of the accretion disk 
temperature distribution more dominant at optical wavelengths. As shown in Table 1 and Figure 3, 
this model predicts essentially no variation of line width with lower level excitation potential and a 
very slight dependence on wavelength in the optical region. (Note that PH08 predict a much larger 
effect of line width on excitation potential than Welty et al. (1992) for what should be essentially 
the same disk model; the reason for the large discrepancy is unclear.) In any case, measurement of 
rotation is best done through cross-correlation using suitable templates, as many of the lines used 
by PH08 are blends and introduce very large scatter into the model predictions (see Table 1 and 
Figure 3). 

PH08 also noted that their disk model predicts very strong TiO absorption bands at ~ 7054 
and 7087 A which are not observed. However, as shown in Figure 2, our lower-temperature disk 
model does not predict strong TiO absorption bands in this region. In addition there is evidence for 
7087 A bandhead absorption in our MIKE spectra, at the level predicted by our disk model. Once 
again this difference in the predicted disk model spectra arises simply by reducing the maximum 
disk temperature, which increases the importance of hot inner disk relative to the outer cool disk 
at the wavelength of observation. 

It has long been known that many optical line profiles in FU Ori are less double-peaked than 
predicted by simple quiescent disk models (Hartmann &; Kenyon 1985). There are, however, alter- 
native possibilities to explain the profiles which do not demand abandonment of the accretion disk 
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hypothesis. If, as currently thought, ionized disks accrete through the action of the magnetorota- 
tional instabihty (MRI; Balbus &: Hawley 1998), such disks must be turbulent. The disk models 
with resolved vertical structure computed by Miller & Stone (2000) predict that turbulence driven 
by the MRI in the central layers of the disk produces waves which propagate outward and shock in 
the upper layers. It would be surprising if the MRI did not produce significant turbulent motions in 
the upper atmospheric layers of the disk, which would tend to wash out the double profile structure. 
Hartmann, Hinkle, & Calvet (2004) found that that some mildly supersonic turbulence was needed 
to explain the ^^CO first-overtone lines of FU Ori. 

It should be noted that the standard steady disk structure may not be completely applicable 
in the innermost disk. Standard thin-disk models predict that accretion onto a slowly rotating 
star should give rise to boundary-layer emission with roughly half the system luminosity; this is 
not observed in FU Ori (Kenyon et al. 1989). Popham et al. (1996) considered disk models which 
suppress boundary layer radiation; such models exhibit less doubled line profiles in inner disk 
regions, largely because the angular velocity of the disk departs from Keplerian values near the 
inner disk boundary. 

To explain the optical spectrum of FU Ori with a central star, the star would be required to 
have essentially the same total system luminosity L ~ 23OL0, and would need to have a radius 
roughly twice the inner radius of the disk model, R ~ IOR0 (Zhu et al. 2007). Assuming Keplerian 
rotation for the infrared disk, Zhu et al. estimated a central mass M ~ 0.3Mq. Such a star 
cannot be an isolated product of stellar evolution, as it has an implausibly short Kelvin-Helmholtz 
contraction time ~ GM^R^^L^^ ~ 1200 yr. The energy to power the star would have to come from 
disk accretion, which would also potentially explain the outburst (Larson 1983). However, as the 
ratio of optical-to-infrared rotational velocities is consistent with a Keplerian profile, this implies 
that any central star would have to be rotating nearly at breakup; this means that the assumption 
of solid-body rotation in the PII08 model is unlikely to be correct. It is also unclear whether the 
accretion of a large amount of hot disk material would add enough angular momentum to spin up 
the outer layers of the star to breakup velocity as it expanded the outer atmosphere. 

In summary, the accretion disk model for FU Ori provides a coherent explanation of the ob- 
served spectral energy distribution and differential rotation over more than a decade in wavelength. 
The slow rotation observed at 5/im supports our previous result that the high mass accretion rate 
disk could extend to 0.5 — 1 AU, which is significantly larger than that predicted by the pure 
thermal instability theory (jBell &: Linlll994l'). On the other hand, the theory incorporating both 
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cessfully predicts AU scale high mass accretion rate inner disk during outbursts (IZhu et al.ll2009l ). 
With the advent of more powerful computers and sophisticated magnetohydrodynamic codes, and 
the assumption of MRI-driven accretion, it should be possible to explore the possibility that atmo- 
spheric turbulence and /or nonstandard inner thin disk structure can explain details of the optical 
line profiles. 
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Table 1: HWHD of lines selected by PH08 as measured from Our synthetic disk spectrum 
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"HWHD means "half- width at half-depth" as illustrated and defined in Figure 3 of PH08. Grade indicates the amount 
of line blending due to rapid rotation: grades and 1 signify badly blended lines, grade 2 identifies less-blended lines 
that are difficult to use, and grade 3 means that the line is moderately blended so that the true HWHD may be 
slightly smaller than our measured value to within 10 km s~^. Grade 4 means the line is slightly blended with an 
errorbar of only 3 km s~^, while grade 5 means that the line is blend-free. 
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Fig. 1. — Observed 4.9 um Phoenix spectrum of FU Ori (solid curve) compared with our synthetic 
disk model spectrum (upper dotted curve). While the wavenumber is shown at the bottom, the 
wavelength in microns is shown at the top of the figure. The lower dotted curve shows the model 
spectrum observed pole-on, which removes the rotational broadening, illustrating the need for rapid 
rotation to match the observations (sec the text). The strong lines are identified in the lower model 
spectrum. The arrow points to the strong CO line which is missing in our model. 
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Fig. 2. — Comparison of the observed MIKE spectrum of FU Ori near 7080 A (solid curve) with 
the synthetic disk spectrum (dotted curve). Compared with the continuum flux ~ 7080 A, there is 
some evidence for weak TiO 7087 A band absorption in both the observed and the model spectrum 
(see text). 
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Fig. 3. — Dependence of the model line half-width at half-depth (HWHD) on excitation potential 
(EP) (left) and the wavelength (right), taken from Table 1. We show only the grade 4 and 5 line 
widths for clarity (see Table 1 for their definitions). 



